Coulomb Blockade in Graphene Nanodisks 
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Graphene nanodisk is a graphene derivative with a closed edge. The trigonal zigzag nanodisk 
with size A'^ has A''-fold degenerated zero-energy states. We investigate electron-electron interaction 
effects in the zero-energy sector. We explicitely derive the direct and exchange interactions, which 
are found to have no SU(A'^) symmetry. Then, regarding a nanodisk as a quantum dot with an 
internal degree of freedom, we analyze the nanodisk-lead system consisting of a nanodisk and two 
leads. Employing the standard Green function method, we reveal novel Coulomb blockade effects in 
the system. The occupation number in the nanodisk exhibits a peculiar series of plateaux and dips, 
reflecting a peculiar structure of energy spectrum of nanodisk without SU(A^) symmetry. Dips are 
argued to emerge due to a Coulomb correlation effect. 



I. INTRODUCTION 

Graphene-related materials has opened a new excit- 
ing field of research in nanoelectronics.— i^i^ In partic- 
ular, graphene nanoribbona^'^'^'^'^ i^^i^^i^^i^^ have at- 
tracted much attention due to a rich variety of band 
gaps, from metals to wide-gap semiconductors. Another 
class of graphene derivatives are graphene nanodisks i^^i^^ 
They are nanometer-scale disk-like materials which 
have closed edges. Graphene nanodisks can be con- 
structed by connecting several benzenes, some of which 
have already been manufactured by soft-landing mass 
spectrometryii2ii2iiiS Nanodisks as well as nanoribbons 
would be promising candidates of future electronic 
nanodevices.— 

There are varieties of nanodisks, among which trig- 
onal zigzag nanodisks are prominent in their electronic 
property because there exist half-filled zero-energy states. 
This novel property was revealed first based on tight- 
binding modeU^ and confirmed subsequently by first- 
principle calculationsji^iii We introduce the size param- 
eter N for trigonal zigzag nanodisks as illustrated in 
FigHUa). Then, there exists TV-fold degenerated zero- 
energy states in the trigonal zigzag nanodisk with size 
N . We have already arguedi^ that spins make a ferro- 
magnetic order and that the relaxation time is quite large 
even if the size N is very small. 

In this paper we make an investigation of electron- 
electron interaction effects in the zero-energy sector of 
the graphene nanodisk consisting of N states. We de- 
rive explicitly the direct and exchange interactions, where 
there is no SU(iV) symmetry. We estimate the spin stiff- 
ness. It is found to be as large as a few hundred meV, 
which means that a nonodisk is indeed a rigid ferromag- 
net. Then we analyze the nanodisk-lead system based on 
the standard Green function method, where a nanodisk is 
connected to right and left leads [Fig[ljb)]. The analysis 
can be carried out quite in the same way as in conven- 
tional quantum dots. The present system turns out to 
be akin to a single dot system with an internal degree of 
freedom. We find a series of Coulomb blockade peaks in 
the conductance as a function of the chemical potential. 
Furthermore, the occupation number in the nanodisk ex- 



hibits a pecuHar series of plateaux and dips, refiecting 
a peculiar structure of energy spectrum without SU(iV) 
symmetry. 



This paper is organized as follows. In Sec. m we sum- 
marize the basic nature of trigonal zigzag nanodisks. In 
Sec. Illli we study direct and exchange Coulomb inter- 
actions in the zero-energy sector. In Sec. IIV^ formulat- 
ing the nanodisk-lead system, we construct the effective 
Hamiltonian near the half filling. In Sec. [Vl we investi- 
gate the Coulomb blockade in the nanodisk-lead system, 
and calculate numerically the occupation number in the 
nanodisk and the Coulomb blockade peaks as a function 
of the chemical potential. Sec. IVjis devoted to conclud- 
ing remarks. 



(a) nanodisk n=o 

N=1 
N=2 
N=3 
N=4 

(b) nanodisk-lead system 





FIG. 1: (a) Geometric configuration of trigonal zigzag nan- 
odisks. It is convenient to introduce the size parameter A*' in 
this way. The 0-trigonal nanodisk consists of a single Ben- 
zene, and so on. The number of carbon atoms is given by 
Nc = A''^ + &N 4- 6. (b) Illustration of the nanodisk-lead sys- 
tem. A nanodisk is connected to the right and left leads by 
tunneling coupling. 
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II. GRAPHENE NANODISKS 

Graphene nanodisks are graphene derivatives which 
have closed edges. The Hamiltonian is defined by 

H = Y,e4c, + Y.^,jc\c,, (2.1) 

where is the site energy, is the transfer energy, and 
c| is the creation operator of the tt electron at the site 
i. The summation is taken over all nearest neighbor- 
ing sites (j, j). Owing to their homogeneous geometrical 
configuration, we may take constant values for these en- 
ergies, £i = Ep and tij = t. We choose t = 3eV as a 
phenomenological parameter—. 

In a previous worki^, we have investigated the elec- 
tronic and magnetic properties of graphene nanodisks 
with various sizes and shapes in quest of zero-energy 
states or equivalently metaUic states. The emergence of 
zero-energy states is surprisingly rare. Among typical 
nanodisks, only trigonal zigzag nanodisks have degener- 
ate zero-energy states and show metallic ferromagnetism, 
where the degeneracy can be controlled arbitrarily by de- 
signing the size. 

It is convenient to introduce the size parameter N for 
trigonal zigzag nanodisks as in Fig[T][a). The size- TV nan- 
odisk has A^-fold degenerated zero-energy states^^, where 
the gap energy is as large as a few eV. Hence it is a good 
approximation to investigate the electron-electron inter- 
action physics only in the zero-energy sector, by project- 
ing the system to the subspace made of those zero-energy 
states. 




FIG. 2: The 3 degenerate zero-energy states of the trigonal 
nanodisk with size A'^ = 3. The solid (open) circle denotes 
that the amplitude Lof is positive (negative). The amplitude is 
proportional to the radius of circle. It is seen that no electrons 
are present in B sites. 

Let \ fa) be the zero-energy state, a = 1, 2, • • • , A^. The 
wave function of the state \fa) is expanded as 

U{x)^Y.^tv^{x), (2.2) 



where ipi {x) is the Wannier function localized at the site 
i, and iuf is the probability density to find an electron 
there. The constraint condition is Ylii^fY = 1- 

Diagonalizing the Hamiltonian l|2.ip , we are able to cal- 
culate explicitly the amplitude wf for zero-energy states 
in the trigonal zigzag nanodisk. All of them are found to 
be real. As an example we show them with size TV = 3 in 
FiglSl where the solid (open) circles denote the amplitude 
Lof is positive (negative). The amplitude is proportional 
to the radius of circle. It is seen that no electrons are 
present on B sites, where the graphene honeycomb lat- 
tice is made of A and B sites. 



III. ELECTRON-ELECTRON INTERACTIONS 

We take two states |/a) and I//3), a 7^ /3, each of which 
can accommodate two electrons with up and down spins 
at most. The two-state system is decomposed into the 
spin singlet xs and the spin triplet xt with the normal- 
ized wave functions, 

/ss(x,x') - -^(/„(x)/^(x')-f/„(x')/Mx))xs, 

(3.1a) 

/st(x,x') = i= (/.(x)/^(x') - /„(x')/Mx))xT. 

(3.1b) 

The Coulomb energies are 

(/ss|i^c|/ss) =(7„;3 + Ja/3, (3.2a) 

(/sT|i?c|/sT) =C/a/3- Ja/3, (3.2b) 

with 

C/a/3 = j d\<fy /:(x)/„(x)F(x - y)/|(y)/My), 

(3.3a) 

Ja/3 = j d'xd^y /:(x)/^(x)V^(x - y)/|(y)/o(y), 

(3.3b) 

where V{x — x') is the Coulomb potential; /*(x)/q,(x) 
is the electron number density in the state \fa), and 
f*{x)f/3{K) is the overlap of the wave functions associated 
with the two states and |/^). They are interpreted 
as the direct and exchange energies. 

Applying the above argument to the many-state sys- 
tem, the effective Hamiltonian is derived as 

= J2 Ucpn (a) n (/3) 

Q>/3 

- ^ E ^"/3[4S(«) • S(/3) + n (a) n (/3)], (3.4) 

a>/3 

where n (a) is the number operator and S(a) is the spin 
operator, 

n{a) = d'l{a)dcr{a), S{a) = ^dl{a)T„a'da'{a), 

(3.5) 
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with dcr(a) the annihilation operator of electron with 
spin a =t,i in the state \fa)' t is the Pauli matrix. 
Note that we have included the on-state Coulomb term 
Uaan (a) n (a) in the effective Hamiltonian p.4p . 

We expand Ua0 and Jap in terms of the Wannier func- 
tions, 

Ua/3 =^wrc^>fwf fd^xd^y 

s •' 

X V'*{x)iP]{x)V{yi - y)(pl{y)(pi{y), (3.6a) 

X ^*{x)^^{y)V{^~y)^l{y)^i{x). {3.6b) 

The dominant contributions come from the on-site 
Coulomb terms with i = j = k = I both for the direct 
and exchange energies. We thus obtain 

Uap^Jap=^UJ2i^?^^)', (3-7) 

i 

with 

U = j d'xd^'v ^*{x)^,{x)V(^-y)^*{y)^,{y). (3.8) 

The nearest neighbor Coulomb interaction vanishes since 
there are no electrons in the nearest neighboring sites: 
See FigH 

We have numerically calculated Uap for nanodisks with 
several size A*". For the case of iV = 2, 

Uii = U22 ^ 0.145C/, U12 = 0.0482C/, 

which gives U11/U12 = 3. For the case of iV = 3, 

Uii =U22 = 0.09 44C/, U12 = 0.0315C/, 
U33 =0.0833C/, Ui3 = U23 = 0.0655[/, 

which gives Un/Uu — 3, U13/U12 — 2. For the case of 
TV = 4, 

Uii =0.0556t/, U12 = 0.0228f/, U13 = Uu = 0.0435f7, 
U22 =0.0768C/, U23 = C/24 = 0.0478C/, 
U33 ^U44. = 0.0934f7, C/34 = 0.0311C/, 

which gives U33/U34 = 3, U13/U12 ~ 2. 

We make some remarkable observations. First, the ex- 
change energy is as large as the direct energy, which is the 
order of a few hundred meV. Thus the spin stiffness Jap is 
quite large, implying that nanodisks are rigid ferromag- 
nets. Furthermore, as we have seen numerically, all Jap 
are the same order of magnitude for any pair of a and /?, 
implying that the SU(A'') symmetry is broken but not so 
strongly in the Hamiltonian (|3.4p . Hence, the zero-energy 
sector is described by the SU(N) Heisenberg model as a 
rough approximation. These facts confirm our previous 
resultflfi' that the relaxation time of the ferromagnetic-like 
spin polarization is quite large even if the size of trigonal 
zigzag nanodisks is very small. 



IV. NANODISK-LEAD SYSTEM 

We investigate the Coulomb blockade of nanodisks in a 
similar way to the case of conventional quantum dots. We 
consider the system comprised of nanodisks with right 
and left connected leads [FigdJJb)]. The Hamiltonian of 
the system is written as 

H = Hu+H^ + Ht:, (4.1) 

where Hb is given by p.4p . and 

H^=J2e (k) {4lct + cllct) , (4.2a) 

ka 

E E i^^i^lMc^) + 4(«)cL) 

kij a 

kty a 

The Hamiltonian TJl describes a noninteracting electron 
gas in the leads with e (k) = h'^k'^/2m, while Ht the 
tunneling interaction between the leads and the nanodisk 
with iL(R) the tunneling coupling constant: We have as- 
sumed that the spin does not fiip in the tunneling process. 
It is convenient to make the transformation 

(l:)^K4S)(l) 

with 

i=^J\tLf + \tR\\ (4.4) 

so that the right and left leads are combined into the 
"even" and "odd" leads. The lead Hamiltonian i?L is 
invariant under above transformation, 

H^ = J2^ (fc) (414. + 4l4.) , (4.5) 

ka 

but the transfer Hamiltonian is considerably simplified, 
- ?E E + 4(")cL) ■ (4.6) 

ka a 

It looks as if the tunneling occurs only between the 
"even" lead and the nanodisk. 

In the case of TV = 1, the system is reduced to the 
simple Anderson model of quantum dot 

i^D -E^ii^T^i, (4.7) 

and hence the physics is well known^-. However, for a 
general case, it is hard to study the electronic properties 
of the system because the interaction strengths Uap take 
complicated values. 

Another simple example is given by the system at the 
half-filling with N > 2. As we have argued, the ground 
state is ferromagnet and the relaxation time is very large 
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even though A'' is small. Based on this fact we simplify 
the Hamiltonian l|3.4p by assuming that the nanodisk is 
a rigid ferromagnet with the spin direction |. Then, only 
electrons with spin t are dynamical within the nanodisk. 
We explore the electric properties of such a system. 

The nanodisk Hamiltonian is described only by the 
electrons with the spin t, 



where 



a>l3 



13 



(4.; 



(4.9) 



Here, UajS stands for the Coulomb energy between the 
added up-spin electron in the state a and the background 
down-spin electron in the state Thus, Sa is the poten- 
tial energy increased when one electron is added to the 
state a. 

(a) double dots 



R 



(b) nanodisk 



tl/«) 



R 



FIG. 3: (Color online) (a) Illustration of the standard double- 
dot system. There is intra-hopping between dots, (b) Il- 
lustration of the nanodisk-lead system, where all down-spin 
electrons are filled up in the nanodisk. An up-spin electron 
tunnels from the left lead to one of the states \fa) and then 
to the right lead. 

We first consider the case where one electron is tun- 
nelled from the lead into the nanodisk. Then it is enough 
to consider only the first term in the Hamiltonian l|4.8p . 
We have numerically calculated Eq for some nanodisks. 
For = 2, 



ei = £2 = 0.193C/. 



For = 3, 



ei = £2 = 0.191[/, 



0.2UU. 



(4.10) 



(4.11) 



For N = A, 

£i = 0.165f7, £2 = 0.195C/, £3 = £4 = 0.216C/. (4.12) 

This Hamiltonian looks similar to that of the iV-dot 
system. However, there exists a crucial difference. On 
one hand, in the ordinary A^-dot system, an electron hops 
from one dot to another dot [FiglSja)]. On the other 



hand, in our nanodisk system, the index a of the Hamil- 
tonian runs over the A^-fold degenerated states and not 
over the sites. According to the Hamiltonian l|4.2bp . an 
electron does not hop from one state to another state 
[FiglUb)]. Hence, it is more appropriate to regard our 
nanodisk as a one-dot system with an internal degree of 
freedom. This fact simplifies the analysis considerably. 



V. COULOMB BLOCKADE IN NANODISK 

We investigate the Coulomb blockade in the nanodisk- 
lead system. It has been argued that the conductance is 
given by the formulck22i. 



G 



2e2 



E 



pLpR 



dnp {e) 



(5.1) 



where A (a, £) is the spectral function of the nanodisk 
system, np (£) is the Fermi distribution function. The 
coupling strength is given by 



np (£) \U 



r for |£| < D 
for £ > £1 ' 



(5.2) 



where p (e) is the density of state in the lead at energy £, 
which is almost a constant within the band limit, |£| < D. 

The spectral function is given by the retarded Green 
function as 

^(a,£) = --ImG^(a,£), (5.3) 

TT 

with 

G^(a,i-i') - -ie{t~t'){{d^{t),di{t')}) . (5.4) 

It exhibits the energy spectrum of the states indexed by 
the quantum number a. At the zero-temperature, the 
formula (15.11) is reduced to 



G 



E^(«,m) 



(5.5) 



with ^ the chemical potential, where the conductance is 
simply proportional to the sum of the spectral density. 
This formula can be interpreted as follows. If there is a 
state at the Fermi energy, an electron can tunnel from the 
lead to the nanodisk via the state by baUistic transport. 
If not, an electron needs to go via high energy states with 
finite gap. 

It is a straightforward task to derive the spectral func- 
tion by employing the standard technique^, as we de- 
scribe in the appendix. Here let us summarize the re- 
sult. The spectral function is simply given by the sum of 
Lorentzians, 



A{a,e) = -Y, 



ra)'+r2 

Mr 



(5.6) 
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It has peaks at the energies Sa and + Uap with broad- 
ening r. The height of the peak at energies is propor- 
tional to 1 — (n/3), and at energies ea + Uajj is proportional 
to {np). The occupation number {np) should be deter- 
mined self-consistently by solving 



(iV-l)(: 



la) = j de np {e)A{a^e) 



(5.7) 



We insert the spectral function l|5.6p into the above equa- 



tion, and obtain the linear equation for 

{N - 1) {ua) 

, ("/3> 



tan 1 £fi _ tan 



Ual3 



r 



(5.8) 



Our task is to solve this equation to determine (ua) ■ 

In order to get an overview of the result, we first 
consider the atomic limit, which is the zero tunneling- 
coupHng Hmit (i — > 0), implying that F ^ 0. The equa- 
tions become very simple and we can obtain analytical 
results. This is indeed a good approximation because T 
is very small. In this limit the spectral function (|5.6p 
becomes 

+ {np) S {e - Ea - Uafj)], (5.9) 
and the self-consistent equation l|5.8p becomes 

{N-l){na)=Y.[{l-{np))e{^l-ea) 

+ {np)e{^JL-ea-Uap% (5.10) 

where 9 {x) is the step function: 6* (x) = for x < and 
9 {x) — 1 for a; > 0. It is easy to see that {ua) — 
if /i < Ea, and that — 1 \i ijl > £a+max^[C/Q^]. 

Though the occupation numbers (ria) are nontrivial in 
the other region, it is straightforward to determine them. 

We have shown A{a,£) and ( function of the 

chemical potential /i in FigHl For the case of iV = 3, 
the first plateau emerges at £1 = £2 < ^ < £3 with 
rti = 712 = 2/3, and the second plateau emerges at £3 < 
/i < £1 + U12 with ni = 712 = 713 = 1/2. For the case of 
N — A, the first plateau emerges at £2 < ^ < £3 = £4 
with 712 = 2/3, the second plateau emerges at very small 
region £3 = £4 < < £2 -f U12 with ni = 712 = = 2/5, 
and the third plateau emerges at £2 + C/12 < /i < £2 + ^^23 
with 7i2 = 4/5, 713 = 3/10, and so on. 

We go on to study the case with a finite value of the 
tunneling coupling. We have calculated numerically the 
occupation number {ua) by solving l|5.8p . and then the 
spectral function A{a,e) by returning it to l|5.6p . We 
show the occupation number and the conductance at the 
zero-temperature in FigHJ 
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FIG. 4: (Color online) The occupation numbers {ria) against 
the chemical potential in the case of = = 0. (a) for 
iV = 3, (b) for TV = 4, (c) for iV = 5, and (d) for = 6. Bold 
red lines denote the doubly degenerated occupancy, while thin 
blue lines denote non-degenerated occupancy. 



Coulomb blockade peaks appear sX ^ — Sa and /i ~ 
£a + Uaf3, where new channels open. There are many 
peaks because the symmetry of the interactions is low, for 
example, Un ^ C/12. There would be only two peaks if 
the SU(A^) symmetry were exact. There are no peaks cor- 
responding to the energy £„ -I- Uap + Ua-y because we have 
neglected the double occupancy in the nanodisk state. 

It is intriguing that the occupancy (ua) shows a pecu- 
liar behavior containing dips as a function of the chemi- 
cal potential ^. Note that, in the usual Coulomb block- 
ade, a monotonic behavior without dips are found in (ria) 
against the chemical potential. The occupancy makes a 
drastic change at the energy of the Coulomb blockade 
peak. 

Let us explain the origin of dips by taking an instance 
of the TV = 3 nanodisk [FiglUa)]. First, an electron en- 
ters the states |/i) or I/2) at /i = £1 = £2, where all other 
up-spin states are empty, as results in the increase of the 
occupancy (711) = (7i2)- Similarly, an electron enters the 
state 1/3) at /I = £3 together with the increase of (713). 
However, this is accompanied with a dip in (711) — (712) 
due to a correlation effect: An electron in the state I/3) 
prefers the absence of electrons in the states |/i) and 
1/2) because of the Coulomb repulsion, thus decreasing 
the occupancy (711) = (712)- Next, an electron enters the 
states 1/1) at — El + C/12 in the presence of an up-spin 
electron in the state I/2), as results in the increase of 
(711). In this way, the occupancy (ua) shows a sudden 
increase at £q and at £„ -I- Uap, and a decrease at ep and 
£/3 + Upj for (3 ^ a. 

It is interesting that the peak at £1 -I- C/13 is tiny in 
FiglDJa). It is interpreted as follows. The peak occurs 
when an electron enters the state |/i) in the setting that 
the state I/3) is already occupied by an up-spin electron. 
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(b) N=4 

















1 f 2/3 


J 




In 





0.14 0.16 0.1S 02 0.22 024 0.26 0.28 



6000 






5000 












4000 








3000 








2000 




1 




1000 













(c) N=5 
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FIG. 5: (Color online) The occupation numbers (ria) (up- 
per figure) and the dimensionless conductance G/(^|- ja_|5rJ^ ) 
(lower figure) against the chemical potential, where it is taken 
that = = 0.001. (a) for = 3, (b) for iV = 4, (c) 
for N = 5, and (d) for A'^ = 6. Bold red lines denote the 
doubly degenerated occupancy, while thin blue lines denote 
non-degenerated occupancy. 



Though such a setting is classically impossible because 
El < S3, it can occur quantum- mechanically. Hence the 
process generates a tiny peak. Similar peaks are observed 
at ei + Ui2 and ei + U13 in an instance of the N ~ A 
nanodisk [FiglUb)]. 

In passing we note that the interaction strengths Sa 
and Uaf3 can be determined experimentally by measuring 
Coulomb blockade peaks. 



symmetry. We have revealed a novel series of dips in the 
occupation number on nanodisk together with Coulomb 
blockade peaks in the conductance as a function of the 
chemical potential. Dips are argued to emerge due to 
a Coulomb correlation effect. It is interesting that the 
interaction strengths Sa and Ua0 can be determined ex- 
perimentally by measuring these peaks. 

We have estimated the spin stiffness in nanodisks at 
the half filling, which is found to be as large as a few 
hundred meV. Thus there exists a strong ferromagnetic 
coupling in the ground state of nanodisk. The relaxation 
time is quite large even if the size is small, and the ground 
state can be regarded as a rigid ferromagnet. Hence only 
electrons with spin opposite to the direction of this fer- 
romagnetic state can go through the nanodisk. We may 
regard the nanodisk-lead system as the spin filter pro- 
vided the electron spin does not fiip during the tunneling 
process between the nanodisk and leads. 

I am very much grateful to P. Kim and N. Nagaosa for 
many fruitful discussions on the subject. The work was in 
part supported by Grants-in-Aid for Scientific Research 
from Ministry of Education, Science, Sports and Culture 
(Nos.070500000466). 



APPENDIX A: GREEN'S FUNCTION 

We analyze Green's function (a, t — t'), which is de- 
fined by (|5.4p . Using the Heisenberg equation of motion, 

ihdtda (t) = [da {t) ,Hn + Hi^+ H^], we find 

{ihdt - Ea) {a, t - t') 
^hS (t -t')+J2 Uc^pG^ (a, t)+Y, {k, t~t'), 

k 

(Al) 

where 



G^ {a,t-t') 


= -i9{t 


F^{k,t-t') 




G^ {a,t~t') 





(A2c) 

Note that we have introduced the notation for the 
correlation between the lead and the nanodisk. 

We again use the Heisenberg equation of motion to 
calculate dtF^ {k,t - t'), 



{ihdt - Sk) F^ (fc, t-t')= i*G^ (a, t-t'). 



(A3) 



VI. CONCLUSIONS 

We have investigated correlation effects in the zero- 
energy sector of graphene nanodisks consisting of N 
states. We have derived explicitly the direct and ex- 
change interactions. It is found that there is no SU(A^) 



Making the Fourier transformations of l|Al|) and l|A3|) and 
combining them, we obtain 

{hu; - - {co)) G^ {a,u;)^h + J2 Uc^G^^ {a, lj) , 

(A4) 
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where 



2r 



■In 



D + huj 



D - hw 



(A5) 

is the self-energy, describing virtual electron tunneling 
between the nanodisk and the lead. The real part gives 
a shift of energy, which is to be used to define the new 
on-state energy 

e'^ = ea+ ReS^. 

However, since the energy shift is the order oi eV / D and 
small, we may neglect it. The imaginary part, which is 
independent of the band width Z?, gives the width of the 
energy peak. 

It is necessary to construct {a,t — t') to obtain 
Green's function (a, t — t') by l|A4p . For this purpose 
we again use the Heisenberg equation of motion, 

ihdtG^ (a, t - t') ^hS {t - t') {{nf3 it) do. (t) , 4 {t')]) 

-~ie{t){{-[H,np {t)d^ (i)],4 it')])- 

(A6) 

The commutator in the last term on the right-hand side 
results in the two terms, 



[H, ripda] = [HD,nfjda] + [Ht, npda] , 

where 

[HD,nf3da] 



(AT) 



[HT,nf3da] 



{Sa + Uafi) nfida - ^ Up-^n^n/^da, 

(A8a) 

n/3 ^ tefc -I- ^ f (cldp - d^pCkj da. 

k k 

(A8b) 



Here we neglect higher-order Green functionsSS,!^. First, 
we neglect the last term in (|A8ap . which corresponds 
to neglect the double occupancy states of initial states. 
Next, we neglect the last term in (jASbp . which corre- 
sponds to neglect spin fiips on the nanodisk during the 
tunneling process. In this way we achieve at 

(ihdt -Sa- Uap) G^ (a, t) = hS (t) {np)+J2 tPf {k, t) , 

(A9) 



where we have introducee a new function, 

ik,t- t') = -i6 [t - t') {{np it) ck (t) , 4 it')}) . 

(AlO) 

Thus the equation of motion for G^ (a, t) produces a new 
function. 

It is necessary to analyze the equation of motion of 
(a, t), which leads to the following commutators. 



[H, npCk] = - EfcCfe + [Ht, Ck] , (Alia) 
[HT.npCk] = - inpck -f '^iWkdp - djjCfcj Ck- (Allb) 



Again we neglect the last term, which corresponds to 
neglect a spin-fiip process^^i^^. In this way we achieve at 



[ihdt - ek) Ff {k,t) ^iG^ {a,t) . 



(A12) 



which no longer produces a new function. 

Now, (jAop and l|Al2p make a closed set of the equa- 
tion of motions, which we can solve easily. Making their 
Fourier transformations, we obtain 

[hw-Ea- Uai3 - (w)) G^ (a, Lj)^n (np) . (A13) 

It follows from and (|A13| that 

1 - {np) 



P^a 



{np) 



hw — Ea — Uap — (to) 



(A14) 



with E = hjj. Hence the spectral function is given by the 
sum of Lorentzians, 



A{a,E) 



TT ^ 

P^a 



(i-(»/3))r 

[E-Eaf+T^ 

{np)T 

(£-£„- Uapf + r2 



(A15) 



which is 115.61) in text. 



^ K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, Y. 

Zhang, S. V. Dubonos, I. V. Grigorieva, and A. A. Firsov, 

Science 306, 666 (2004). 
^ K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, M. 

I. Katsnelson, I. V. Grigorieva, S. V. Dubonos, and A. A. 

Firsov, Nature 438, 197 (2005). 
^ Y. Zhang, Y. -W Tan, H. L. Stormer, and P. Kim, Nature 



438, 201 (2005). 
* M. Fujita, K. Wakabayashi, K. Nakada, and K. Kusakabe, 

J. Phys. Soc. Jpn. 65, 1920 (1996). 
^ M. Ezawa, Phys. Rev. B, 73, 045432 (2006). 
^ L. Brey, and H. A. Fertig, Phys. Rev. B, 73, 235411 (2006). 
^ F. Munoz-Rojas, D. Jacob, J. Fernandez-Rossier, and J. J. 

Palacios, Phys. Rev. B, 74, 195417 (2006). 



8 



® Y. -W Son, M. L. Cohen, and S. G. Louie, Phys. Rev. 

Lett., 97, 216803 (2006). 
® V. Barone, O. Hod, and G. E. Scuseria, Nano Lett., 6, 

2748 (2006). 

^° M. Y. Han, B. Oezyilmaz, Y. Zhang, and P. Kim, Phys. 

Rev. Lett., 98, 206805 (2007). 
" Z. Chen, Y. -M. Lin, M. J. Rooks, and P. Avouris, 

|cond-mat/0701599| 
1^ Z. Xu and Q. -S. Zheng, AppL Phys. Lett. 90, 223115 

(2007). 

B. Ozyilmaz, P. Jarillo-Herrero, D. Efetov and P. Kim, 
cond-mat/arXiv:0709.1731. 
" M. Ezawa, Physica Status Solidi (c) 4, No. 2, 489 (2007). 
M. Ezawa, cond-mat/arXiv:0707.0349 (to be published in 
PRB). 

^® J. Fernandez-Rossier, and J. J. Palacios, Phys. Rev. Lett. 
99, 177204 (2007). 



O. Hod, V. Barone, and G. E. Scuseria, cond- 
mat/arXiv:0709.0938. 

C. Berger, Z. Song, X. Li, X. Wu, N. Brown, C. Naud, D. 
Mayou, T. Li, J. Hass, A.N. Marchenkov, E. H. Conrad, 
P.N. First and W. A. de Heer, Science 312, 119 (2006). 
H. J. Rader, A. Rouhanipour, A. M. Talarico, V. Palermo, 
P. Samori, and K. Miillen, Nature materials 5, 276 (2006). 
R. Saito, G. Dresselhaus, and M. S. Dresselhaus, Physical 
Properties of Carbon Nanotubes, Imperial College Press, 
1998, London. 

21 P. W. Anderson, Phys. Rev, 124, 41 (1961). 

22 Y. Meir and N. S. Wingreen, Phys. Rev. Lett., 68, 2512 
(1992). 

23 L. Lacroix, J. Phys. F, 11, 2389 (1981): L. Lacroix, J. 
Appl. Phys., 53, 2131 (1982). 

2* G. Czycholl, Phys. Rev. B, 31, 2867 (1985). 



